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ABSTRACT 

We present a dynamical study of NGC 4697, an almost edge-on, intermediate- 
luminosity, E4 elliptical galaxy, combining new surface brightness photometry, new 
as well as published long-slit absorption line kinematic data, and published plane- 
tary nebulae (PNe) velocity data. The combined kinematic data set extends out to 
~ 5' ~ 4.5i?e and allows us to probe the galaxy's outer halo. 

For the first time, we model such a dataset with the new and flexible x^-made- 
to-measure particle code NMAGIC. We extend NMAGIC to include seeing effects, 
introduce an efficient scheme to estimate the mass-to-light ratio, and incorporate a 
maximum likelihood technique to account for discrete velocity measurements. 

For modelling the PNe kinematics we use line-of-sight velocities and velocity dis- 
persions computed on two different spatial grids, and we also use the individual velocity 
measurements with the likelihood method, in order to make sure that our results are 
not biased by the way we treat the PNe measurements. 

We generate axisymmetric self-consistent models as well models including vari- 
ous dark matter halos. These models fit all the mean velocity and velocity dispersion 
data with /N < 1, both in the case with only luminous matter and in potentials 
including quite massive halos. The likelihood analysis together with the velocity his- 
tograms suggest that models with low density halos such that the circular velocity 
Vc < 200km s~^ at 5i?e are not consistent with the data. A range of massive halos 
with Vc — 250km s~^ at 5i?e fit the PN data best. To derive stronger results would 
require PN velocities at even larger radii. The best fitting models are slightly radially 
anisotropic; the anisotropy parameter (3 =~ 0.3 at the center, increasing to ~ 0.5 
at radii > 2i?e. 

Key words: galaxies: kinematics and dynamics - galaxies: elliptical and lenticular, 
cD - galaxies: individual: NGC 4697 ~ galaxies: halos - cosmology: dark matter - 
methods: iV-body simulation 



1 INTRODUCTION 

The presence of dark matter (DM) has long been inferred 
around spiral galaxies from their flat rotation curves, and 
galaxies are now generally believed to be surrounded by ex- 
tended dark matter halos. Indeed, in the current A-cold dark 
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matter (ACDM) cosmology, galaxies form within the poten- 
tial wells of their halos. The standard picture for the forma- 
tion of elliptical galaxies is through mergers of smaller units. 
Ellipticals should thus also be surrounded by dark matter 
halos. Their halos are particularly interesting because ellip- 
ticals are among the oldest galaxies and are found in the 
densest environments. 

Unfortunately, mass measurements in elliptical galax- 
ies have been difficult because of the lack of a suit- 
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able ubiquitous tracer such as neutral hydrogen rotation 
curves in spirals. In giant ellipticals, there is evidence 
for dark matter from X-ray emission (e.g. | Awaki et alJ 



1 19941 : iLoewenstein fc Whitel [l999l : iHumphrev et all 
and gravitational lensing [e.g. 'Griffiths et al. ' 
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iTreu fc Koopmansll200i : iRusin fc Kochanekii2005D . In more 
ordinary ellipticals, mass estimates come from stellar dy- 
namical studies, which have been limited by the faint- 
ness of the galaxies' outer surface brightness to radii less 
than two effective radi i from the centre, R < 2Re {e.g. 
iKronawitter et al.|[200ol : iThomas et al.ll2007h . These stud- 
ies suggest that the d ark matter contributes ~ 10 — 40% of 
the m ass within Re. l|Gerhard et al.l 1200 ll : ICappellari et al.] 
I2OO6I ). consistent with the lensing results. 

The strong emission line at [OIII]A5007 from plane- 
tary nebulae offers a promising tool to overcome this lim- 
itatio n and to extend stellar kinematic studies to larger 
radii (IHui et al ] | 19951 : iTremblav et al ] |l995l : lArnaboldi et al.l 
Il996l . ll998l )." Also, in the less massive, X-ray faint ellipticals, 
PNe may be the primary tool for constraining the dynamics 
at large radii. Once the PNe are identified, their line-of-sight 
velocities can be obtained from the Doppler shift of the nar- 
row emission line. Interestingly, the derived PNe dispersion 



profi les in the elliptical galaxies NGC 4697 llMendez et al 



200l|) and NGC 821, 33 79 and 4494 (|Romanowskv et al 



2003 : lDouglas et al.|[2007h were found to decline significantly 
with radius outside IRe. Their spherically sy mmetric dy- 
namical analysis led iRomanowskv et al.l (|2003l ) to the con- 
clusion that the se galax i es lac k massive dark matter ha- 
los; however, De kel et al.l (|2005l ) argued that the well known 
mass-anisotropy degeneracy allows for declining dispersion 
profiles even when a standard DM halo is present. 

In the present paper we focus on NGC 4697, a normal 
and almost edge-on E4 galaxy located along the Virgo south- 
ern extension. [Mendcz ot al. (2001) obtained a planetary 
nebula luminosity function (PNLF) distance of 1 0.5 ± 1 Mpc 
from mag nitudes of 531 PNe, and iTonrv et al.1 |200l) mea- 
sured a surface brightness fluctuation (SBF) distance of 
11.7 ± 0.1 Mpc. This fairly isolated galaxy has a total B 
magnitude Bt ~ 10.14 and harbors a central super massive 
black hole (SMBH) of mass 1.2 x lO^Mg (|Pinknev et all 
I2OOOI ). A Sersic law with Re — 66 arcsec gives a good fit to 
the surface brightness profile out to about 4 a rcmin (see Sec - 
tion 2). Based on the disky isophote shapes ICarteil (|l987l ) 
and lGoudfrooii et al. I (I1994') inferred a stellar disk along the 
major axis. The contribution of the disk kinematics to the 
major axis line-of-sight velocity distributions was estimated 
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bv IScorza et all (IT998; '). X-ray observations with ROSAT 



l|Sansom et al.l I2OO0I ) show a lack of large scale hot gas in 
the halo of this g alaxy. Using more recent Chandra data, 
llrwin et al.l (|2000l ) could resolve most of this emission into 
non-uniformly distributed low mass X-ray binary (LMXB) 
point sources, suggesting that NGC 4697 is mostly devoid 
of interstellar gas and perhaps does not have substantial 
amounts of DM. 

Dynamical models of NG C 4697 have bee n cons tructed 
by iBinnev et al.l (|l990l ) and iDeionghe et al. I (|l996l ). both 
based on photometry and absorption line kinematic data 
within ~ IRe- The data were consistent with a constant 
mass-to-light-ratio an d none of these mod els showed evi- 
dence for dark matter. [Mendez et al.l (|200ll ) obtained veloc- 
ities for 531 PNe and derived a velocity dispersion profile out 



to approx imately 4.5Re. Assuming an isotropic velocity dis- 
tribution, iMendez et al.1 (120011 ) found that the PNe velocity 
dispersion profile is consistent with no DM inside 4.5i?e, 
but that DM can be present if the v elocity distr i bution 
is anisotropic. This was also arg ued bv iDekel et al.1 (120051) 
to be the main cause of the finding by IRomanowskv et al.l 
(|2003l ). that their three intermediate lumi nosity galaxies 
lacked significant dark matter halos (but see lDouglas et aU 
120071) . Contrary to these three galaxies, which are nearly 
round on the sky, NGC 4697 is strongly flattened and likely 
to be nearly edge-on, thus easier to model since shape de- 
generacies are much less severe. 

In the light of this, it is important enough to perform a 
further careful analysis of this galaxy. In this paper we con- 
struct dynamical models of NGC 4697 with the very flexible 
NMAGIC particle code , making use o f new and published 
slit kinematics and the iMendez etHI (|200ll ) PN data. The 
paper is organized as follows. In Section [5] we describe our 
new observational data and all other data that are used 
for the dynamical models. In Section [3] we give a brief ex- 
planation of the NMAGIC modeling technique. We extend 
the method to include seeing effects, introduce an efficient 
scheme to estimate the mass-to-light ratio, and show how 
discrete velocity measurements may be taken into account. 
In Section |4] we construct isotropic rotator models to test 
and calibrate the method, preparing for the dynamical mod- 
eling of NGC 4697 which is then performed in Section [S] 
Finally, the paper closes with our conclusions in Section [S] 



2 OBSERVATIONAL DATA 

In this section, we describe the observational data used in 
the present study, including new CCD photometry and new 
long-slit absorption line kinematics. We also describe here 
the procedure employed for the deprojection of the photo- 
metric data. I n the following we a dopt a distance 10.5 Mpc 
to NGC 4697 l|Mendez et al.ll200ll ). 



2.1 Surface Photometry 

The R-band data used in the present work were taken in 
AprU 2000 as part of the ESQ Proposal 064.N-0030 (P.I. R.P. 
Saglia) at the Wide Field Imager on the ESO-MPIA 2.2m 
telescope. Six 5 minutes, dithered exposures where taken in 
sub-arcsec seeing conditions. After the usual data reduction 
procedures (performed using the IRAF task mscred), the 
data were tabulated as radial profiles of surface brightness fi, 
ellipticity e, position angle PA and Fourier shape coefficients 
(|Bender fc Moellenhofil Il987l ). The surface brightness was 
calibrated using the R band photometry of iPeletier et al.l 
h99(t ). Systematic errors due to sky subtraction {f-Lsky — 
20.18 Rmag arcsec"^) are always smaller than 10%. Fig. 
lAll shows the derived profiles and Table |X] gives them in 
tabular form. The isophotes of NGC 4697 do not show any 
appreciable twist in PA and have a positive 04 coefficient in 
the galaxy's inner parts , which is wcU-cxplaincd by a near- 
edge-on embedded disk (jScorza fc B ender 1 995). The galax y 
has some dust in the inner regions (jPinknev et al.ll2003l ). 
but the R-band observations are relatively unaffected by it. 
The outer isophotes are progressively slightly off-centered. 
A Sersic fit to the surface brightness profile results in Sersic 
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Figure 1. Comparison of tlie pliotometry of NGC 4697 with 
reprojected tliree-dimensional luminosity models. The data points 
correspond to the observed photometry for NGC 4697 (for the 
deprojection, the and ag values have been set to constants 
beyond ~ 3i?e)- The solid line shows the edge-on deprojected 
model reprojected, the dashed line the i = 80° deprojected model 
reprojected, and the dash-dotted line the i = 67° model. The 
panels show, from top to bottom, surface brightness nn, ellipticity 
e and the isophotal shape parameters a4 and ag. 



index n = 3.5 and effectiv e radius Re = 66 a rcsec. The 
older value of 95 arcsec from lBinnev et all was based 

on photometry reaching only 120 arcsec; thus we use Re = 
66 arcsec in the following. For a distance of 10.5 Mpc this 
corresponds to 3.36 kpc. 



2.2 Deprojection 

For our dynamical study, we will fit particle models to 
the deprojected luminosity density using NMAGIC, cf. Sec- 
tion [3l To obtain the luminosity density we need to de- 
project the surface brightness. This inversion problem is 
unique only for spherical o r edge-on axisymmetric systems 
iBinnev fc Tremaind [l987l ). For axisymmetr ic systems in - 
clined at an angle i, the Fourier slice theorem (lRvbickilll987l ) 
shows that one can recover information about the density 
only outside a "cone of ignorance" of opening angle 90° — i. 
The deprojection of photometry for galaxies with i signif- 
icantly less than 90° can thus be signi ficantly in error be- 
cause of u ndetermined konus densities (|Gerhard fc BinnevI 
ll996l ; [Koc hanek fc Rybicki 1996). 

Fortunately, NGC 4697 is seen almo st edge-on and 
hence does not suffer from this ill-condition. iDeionghe et al.l 
1 19961 ) observed a nuclear dust lane with a ring-like appear- 
ance, elongated along the major axis of NGC 4697. Assum- 
ing that the ring is settled in the equatorial plane, they es- 
timated an inclination angle i = 78° ± 5° . Applying a disk- 
bulge decomposition to th e isophote shapes of th e galaxy, 
and assuming a thin disk, IScorza fc Bendeil I^1995^ derived 
an inclination i — 67°. This was updated by Scorza et al.l 



Figure 2. Isophote parameter profiles for the NGC 4697 photom- 
etry compared with those of the deprojected luminosity models 
shown in Figure [Tl but now as seen in edge-on projection. NGC 
4697 is coded as the solid line, whereas the i = 80° and i = 67° lu- 
minosity models are presented by the dashed and the dash-dotted 
lines, respectively. 



(|l998h to i = 70°. IScorza fc Bendeij (|l995l ) also estimated 
the velocity dispersion aa of the disk component from the 
major axis line-of-sight velocity distributions. From their 
plots we estimate ad — 95km at the half-mass radius of 
the disk, rn — 13". Assuming that the vertical dispersion in 
the disk is comparable, we can estimate the intrinsic flatten- 
ing of the disk oc u^/Ve — 0.2, using the measured rotation 
velocity. A disk with intrinsic thickness h/R ~ 0.2 would 
give the same isophote distortions for inclination i ~ 80° as 
a thin di sk with ^ = 67 — 70°, in agreement with the value 
found bv lDejonghe et al.l 1 19961 ). 

We have deprojected NGC 4697 for inclinations i — 
90°, i = 80° and i = 67°, using the method and pro- 
gram of iMagorriaiil (| 19991 ). This algorithm uses a maxi- 
mum penalized likelihood method with a simulated anneal- 
ing scheme to recover a smooth three-dimensional luminos- 
ity density distribution which, when projected onto the sky- 
plane, has minimal deviations from the observed SB. The 
three-dimensional luminosity density so obtained extends 
beyond the radial range of the data, where the penalized 
likelihood scheme favours an outer power-law density pro- 
file. Figure [T] compares the observed photometry with the 
three deprojected luminosity models reprojected on the sky. 
In the range from 0.2Re to 2Re the i = 67° deprojection 
yields a significantly less good fit to the observed 04 and 
ae. Figure [2] compares the radial run of the isophotal shape 
parameters for the i — 80° and i = 67° luminosity models 
projected edge-on, with the observed photometry of NGC 
4697. The i = 80° deprojection produces again somewhat 
better results. It is also more physical because it allows for a 
finite thickness of the inner disk of NGC 4697, as discussed 
above. Hence we will adopt it for the dynamical study of 
NGC 4697 to follow. 
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2.3 Kinematic Data 

2.3.1 Stellar-absorption line data 

Long slit absorption line kinematics within ^ R e have 
bee n reported, amongst oth er wo rks, by Binnev et alJ (|l990l ) 
and iDeionghe et al.l (|l996h . The lBinnev et all (jl99d ') kine- 
matic data (BDI data) consist of line-of-sight velocity v and 
velocity dispersion a profiles along the major axis, along 
slits 10" and 20" parallel to the major axis, along the mi- 
nor axis, and along a slit 22" parallel to the minor axis. 
They are derived u sing the Fourier Quot ient (FQ) m ethod 
ijllhngworth fc Schec hter 1982). Deionghc ^eTal] ljl996l l have 
published further v and a measurements (DDVZ data) at 
various position angles, also measured with the FQ method. 

Along the major and minor axes we have derived ad- 
ditional line-of-sight velocity distribution (LOSVD) kine- 
matics from the high S/N inte g rated absorption line spec- 
tra obtained by iMendez et"all (120051 ) with FORS2 at the 
VLT, a slit width of 1 arcsec and seeing 1" — l."5. We 
refer to this paper for a description of the data acqui- 
sition and reduction. The LOSVDs were measured us- 
in g the Fourier Corre lation Quotient (FCQ) m ethod, as 
in EenderetiD (|l994l ) and iMehlert et all (|2000l ). and the 
K3III star HD132345 as a template. From these LOSVDs, 
profiles of V, a, fea and hj, the Gauss- Hermite coeffi- 
cients ijGerhardl [l993l : | van der Marel fc Franxl [l993l l . were 
obtained; these are shown in Figure [S] Tables lAll and IA2I 
give the data in tabular form. The statistical errors derived 
from Monte Carlo simulations are minute and much smaller 
than the rms scatter observed between the two sides of the 
galaxy. These differences are particularly obvious along the 
major axis in the r adial range 10-20 arcsec. As noted in 
iMendez et all l|2005l ). in this region we detect patchy [GUI] 
emission that is affecting to some extent the kinematics. 
Judging from the asymmetries in the kinematics on both 
sides of the galaxy, we estimate the residual systematic er- 
rors affecting the data, which amount to ~ 3 km/s in V, 
~ 3.5 km/s in a, ~ 0.02 in /13 and /14. 

In the following we discuss the comparison between the 
kinematic data derived here and the kinematics published in 
the literature. Note that part of the differences seen below 
arise from the different observational setups (along the ma- 
jor ajcis different slit widths probe different relative amounts 
of the central disk structure present in the galaxy and re- 
lated [QUI] emission regions) and methods used. In partic- 
ular, the FQ method fits Gaussian profiles to the LOSVDs, 
ignoring the higher order Gauss-Hermite terms. In general 
the systematic effect on the mea sured mean velocity and ve - 
locity dispersion profiles is small (|van der Marel et al.lll99^ . 
When applied to our major axis data set, FQ gives system- 
atically slightly smaller mean velocities and larger velocity 
dispersions. 

In the inner 10 arcsec along the major axis we confirm 
th e clear kinematic si gnatu re of the central disk discussed 
by IScorza &: Bendeil (|l995l . SB) and agree well with their 
mean velocities and velocity dispersions derived also with 
the FCQ method and 1.8 arcsec slit width (see Figure 31 
left, cyan points). Along the same axi s we find good agree- 
ment with iFried fc Illingwor"thl (|l994l ) (FI, yellow, derived 
using the FQ method and 1.1 arcsec slit width). Overall, 
the BDI data (Figure [4] red points, 2.6 arcsec slit width), 
agree well with our data, although at 15" along the ma- 
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Figure 5. Schematic view of the slit setup used to construct the 
particle models. BDI slits are shown in black and our slits are 
shown in red. The ellipse has a semi-major axis of length Rg and 
axis ratio q = 0.6. 



jor axis, the two data sets differ systematically. The a pro- 
files of DDVZ (blue points, 0.7 arcsec slit width) differ sig- 
nificantly in the sense that at small semi-major axis dis- 
tances the DDVZ a is increasing with radius but our a is 
decreasing. Finally, the right panel of Figure [4] compares 
our data along t he minor axis with th e data sets of BDI 
(red points) and iKoprolin et al.l (|2000l ) (KZ, green points, 
2 arcsec slit width), who use the Fourier Fitting method 
of Ivan der Marel fc Franxl ([l993). Both agree within their 
respective (larger) errorbars. 

Based on the radial extent and quality of the different 
datasets, and taking into account the discussion above, we 
have decided to use in the subsequent modeling only our 
data combined with BDI. For our kinematic data, the errors 
in V and a are < 0.5km s~^, which is small compared to 
the scatter in the data. This suggests, as already discussed 
above, that systematic errors dominate. For the modelling 
we have therefore replaced these errors with the smallest 
errors in v and a of the BDI data along the major and minor 
axis, respectively (5 — 7km s^^). Similar arguments hold for 
hs and /i4 and we set their errors to 0.01. In addition, the 
coefficients along the minor axis scatter significantly around 
zero while the minor axis velocities are consistent with zero; 
thus we replace these values by /13 — 0.0. 

Figure [S] gives a schematic view of the arrangement of 
the kinematic slits used in the modelling process. 



2.3.2 Planetary nebula velocities 

Planetary nebulae (PNe) are dying stars that emit most of 
their light in a few narrow lines of which the [OIII]A5007 
is the most prominent one. The PN population in elliptical 
galaxies is expected to arise from the underlying galactic 
population of old stars and hence the PNe can be used as 
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Figure 3. The kinematics along the major (left panel) and minor (right panel) axis of NGC 4697. The filled and starred symbols refer 
to the data folded along the axes. 
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Figure 4. Comparison between the different absorption line kinematics along the major (left panel) and the minor (right panel) axis. 
Black: our data; red: BDf; blue: DDVZ; cyan: SB; green: KZ; yellow: Ff. 
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3 NMAGIC MODELS 

In this section we give a brief introduction to NMAGIC 
and present a few exte n sions to the method de scribed in 
Ide Lorenzi et alj l|200i1 ). ISver fc Tremaind (|l99i) invented 
a particle-based method for constructing models of stellar 
systems. This "made-to-measure" (M2M) method works by 
adjusting individual particle weights as the model evolves, 
until the N-particle system rep roduces a set of target con- 
straints, ' de Lorenzi et al.] (|2007l ) improved the algorithm to 
account for observational errors and to assess the quality of a 
model for a set of target data directly, using the standard 
statistics in the function to be maximized upon convergence 
of the weig hts (x^M2M). NMAGIC is a parallel implemen- 
tation of the improved x^M2M algorithm. 



Figure 6. The positions and velocities of the cleaned PN sample 
of 381 PNe. The lines indicate the different grids used for binning 
the PNe. Details are given in the text. 



kinem atic tracers for the stellar distribution. iMendez et al.l 
l|200ll ) detected 535 PNe in NGC 4697 and were able to 
measure radial velocities for 531 of these with a typical error 
of 4 km s~^ 

ISambhus et al.l (|2006l ) analyzed the correlations be- 
tween the magnitudes, velocities and positions of these 531 
PNe and found kinematic evidence for more than one PN 
sub-population in NGC 4697. In addition to the main PN 
population, they found evidence for a population of prefer- 
ably bright PNe which appeared to be not in dynamical 
equilibrium in the galactic potential. To remove these pos- 
sible kinematic co ntaminants, and to also ensure complete- 
ness for R > Re ('Mendez et al.ll200ll '). we discard all PNe 
with magnitudes outside the range 26.2 < m(5007) ^ 27.6. 
The positions and velocities of the remaining 381 PNe are 
shown in Figure [6] In the following, we use a doubled sam- 
ple of 762 PNe for our analysis, obtained by applying a 
point-symmetry reflection. Every PN with observed posi- 
tion coordinates {x, y) on the sky and line-of-sight velocity 
vpM is reflected with respect to the center of the galaxy to 
generate a new PN with coordinates {—x, —y, —vpn)- Such 
point-symmetric velocity fields are expected for axisymmet- 
ric and non-rotating triaxial potentials. Moreover, this re- 
flection will help to further reduce any PN sub-population 
biases which might still be present. 

We compute v and a on two slightly different spatial 
grids, subtracting 40 km s~^ in quadrature from all PN ve- 
locit y dispersions to acco unt for the measurement uncertain- 
ties l|Mendez et al.|[200ll ). We use the spatial bins defined by 
the solid lines displayed in Figure[6]to obtain data set PNDl, 
which is shown together with the models in Section [S] The 
second data set, PND2, is computed using the same grid 
but replacing the outermost ellipse by the dashed ellipse 
with semi-major axis a = 280". This second grid is used to 
make sure that the dynamical models we generate are not 
affected by the way we define the outermost bins. 



3.1 Luminous and dark matter distribution 

3.1.1 Luminous mass 

We assume that the luminous mass distribution of NGC 
4697 follows the deprojected luminosity density. The mass 
density of the luminous matter is then given by = Tj, 
with mass-to-light ratio T and luminosity density j repre- 
sented by the discrete ensemble of particles with positions 
Xi and luminosities Li. 



3.1.2 Dark halo potential 

The prevailing cosmological paradigm predicts that galax- 
ies have massive, extended dark matter halos. Numeri- 
cal cold dark matter (CDM) simulations reveal univer- 
sal hal o density profi l es wi t h steep central den sity cusps 
[e.g. iNavarro et al] 1 19961 : iMoore et all Il999l l. On the 
other hand, observations of many dwarf and low-surface 
brightness galaxies find shallower inner d ensity cores {e.g. 
Ide Blok et al]|2003l : iMcGaugh et"al]|2007l ). Here our aim is 
not to determine the detailed shape of the dark matter halo 
in NGC 4697, but rather to first see whether the PN veloci- 
ties allow or require any dark matter at all in this galaxy. To 
answer this question we will investigate a one-dimensional 
sequence of potentials whose circular velocity curves vary 
at large radii between the near-Keplerian decline expected 
when the mass in stars dominates, and the nearly flat shapes 
generated by massive dark halos. Thus for our dynamical 
studies of NGC 4697, we represent the dark matter halo by 
the logarithmic potential (|Binnev fc Tremain&,1987. ) 



which is generated by the density distribution 
Pd{R', z 



(1) 



{2ql + l)4 + R'^ + 2{l-^q-')z'' 



AivGq; 



(2) 



where vo and ro are constants, is the flattening of the 
potential, and R' and z' are cylindrical coordinates with 
respect to the halo's equatorial plane. When q^ < l/\/2 
the density becomes negative along the z' axis. The density 
given in equation [2] has a shallow inner density profile, but 
since we are mainly interested in the circular velocity curve 
in the outer halo of NGC 4697, this is inconsequential: it 
is always possible to reduce the stellar mass-to-light ratio 
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in exchange for an additional centrally concentrated dark 
matter cusp. 



3.1.3 The total gravitational potential 

The total gravitational potential is generated by the com- 
bined luminous mass and dark matter distributions and is 
given by 



(3) 



where <j)* is generated by the A'^-particle system assuming a 
constant mass-to-light ratio for each stellar partic le. We es- 
timate (|>^, via a spherical ha rmonic decomposition jSellwoodl 
l2003l : lde Lorenzi et al.ll2007l ). The stellar potential is allowed 
to change during a NMAGIC modelling run, but the dark 
matter potential is constant in time and is given by equa- 
tion lU. The particles are integrated in the global potential 
using a drift-kick-drift form of the leapfrog scheme with a 
fixed time step. 



3.2 Model observables 

Typical model observables are surface or volume densities 
and line-of-sight kinematics. An observable j/j of a particle 
model is computed via 



N 



WiKj [zi{t)] , 



(4) 



where Wi are the particle weights, Zi are the phase-space 
coordinates of the particles, i = 1, ■ ■ • , A^, and Kj[zi{t)] is a 
kernel corresponding to yj. We use units such that the lu- 
minosity Li of a stellar particle can be written as Li — Lwi, 
where L is the total luminosity of the model galaxy. We 
use temporally smoothed observables to i ncrease the effec- 
tive n u mber of particles i n the system, c/. ISver fc Tremain^ 
(| 19961 ) : lde Lorenzi et all (|2007l l. 



3.2.1 Luminosity constraints 

For modeling the luminosity distribution of NGC 4697 one 
can use as observables the surface density or space density 
on various grids, or some functional representations of these 
densities. We have chosen to model a spherical harmonics 
expansion of the deprojected luminosity density. We deter- 
mine the expansion coefficients Aim for the target galaxy 
on a 1-D radial mesh of radii rt. The spherical harmonic 
coefficients for the particle model are computed via 



ai,n,k = L'y^^jki'^Yi"^ {6i, tfii)wi 



some constraints on which observables can be used. For 
kinematics, suitable observables are the luminosity-weighted 
Gauss-Hermite coefficients or the luminosity-weighted veloc- 
ity moments. We implement them as follows. 

Spectroscopic data The shape of the line-of-sight veloc- 
ity distribution (LOSVD) can be expressed as a truncated 
Gauss-Hermite series and is then characterized by V, a and 
h„ (n > 2), where V and a are free parameters. If V and 
(T are equal to the parameters of t he best fitting Gaussian 
to th e LOSVD, t hen h i = /i2 = (|van der Marel fc Fra"iu3 
| l993j ; iRix et al.l Il997h . The luminosity- weighted Gauss- 
Hermite coefficients are computed as 



bn,p = Ip hn.p = 2y/TlL SpiUn{l^pi)wi, 
i 

with 

l^pi — {Vz,i — Vp) jOp. 



(6) 



(7) 



Here Vz,i denotes the line-of-sight velocity of particle i, Ip 
is the luminosity in cell Cp, Vp and Op are the best-fitting 
Gaussian parameters of the target LOSVD in cell Cp, and th e 
dimensionless Gauss-Hermite functions are (|Gerhardlll993 ) 

u^(y) = (2"+ Vn!) exp (-i^Vs) • (8) 

Hn are the standard Hermite polynomials and Spi is a selec- 
tion function which is one if particle i is in cell Cp and zero 
otherwise. The errors in hi and h2 can be computed from 
those of V and a via 

AH.^-j=^ (9) 
V2 o" 

and 



Aho = 



1 Act 



(10) 



valid to first order (Ivan der Marel fc Franx|[l993 : [r ix et al.l 
Il997f ). Since we use the observed Vp and ap from a Gauss- 
Hermite fit to the LOSVD as expansion parameters for the 
model line profiles, the final fitted hi and /12 of a model will 
be small, and so we can also use relations ((9} and (|10|) to 
compute the model V and a from Vp and Op. 

Spatially bmned PNe data We have computed mean PN 
velocities and velocity dispersions for the ellipse sector bins 
shown in Figure |6] The ellipticity of the grid corresponds to 
the mean ellipticity of the photometry. As suitable observ- 
ables we take the luminosity-weighted velocity moments in 
these bins, which are computed as 



(5) 



(11) 



where L is the total luminosity of the model galaxy, Wi 
the particle weights, Yj™ the spherical harmonic functions 
and yki"^' is a selection function, which maps the parti- 
cles onto the radial m esh using a cloud-in-cell scheme (see 
Ide Lorenzi et al.ll2007l '). 



where Vz,i is the velocity along the line-of-sight of particle i 
and Spi is a selection function, which is equal to one if par- 
ticle i belongs to the bin segment under consideration and 
zero otherwise. In the following, we use only the moments 
v„ and v„. 



3.2.2 Kinematic constraints 

Since in the x^M2M algorithm the kernel in equation Q 
cannot depend on the particle weights themselves, this puts 



3.3 Seeing effects 

To account f or seeing effects we ap ply a Monte Carlo ap- 
proach (e.g. ICappellari et al.]|2006l ) instead of convolving 
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the observables with the PSF. As long as the particles move 
along their orbits no PSF effects need to be taken into ac- 
count, only when the observables of the system are com- 
puted, the effects of seeing may matter. 

When computing an observable including PSF effects, 
we replace the "original" particle at position {xi,yi) on the 
sky plane temporarily by a cloud of Npp pseudo particles. 
The position of a pseudo particle is obtained by randomly 
perturbing {xi,yi) with probability given by the PSF. Note 
that neither extra storage is needed nor additional time 
to integrate the particles along their orbits. Usually, only 
a small number of pseudo particles are needed to model 
PSF effects, even one is often sufficient. This procedure is 
implemented in the kernel Kij as defined in equation Q. 
Th e same kernel th e n ent ers the force-of-change equation, 
cf. Ide Lorenzi et al.1 (|2007| ). 

To test how well PSF effects are modeled using only a 
few pseudo particles we computed mock observations for a 
spherical isotropic galaxy of mass M = 10^" M©, located at a 
dista nce 10 Mpc . The intrinsic density of the galaxy is given 
by a iHernauistI (|l990l ) profile with scale length a = 55.1". 
We assumed a major axis slit of width 2" and a Gaussian 
PSF with FWHM = 4". We computed the LOSVD of the 
target gal axy along the major axis u sing higher order Jeans 
moments l|Magorrian fc Binnevlll994l ) , and compared it with 
the LOSVD of a particle realization of the mock galaxy, ap- 
plying the above procedure. The particle realization of the 
Hernquist model w as generated from an isotrop ic distribu- 
tion function {e.g. iDebattista fc Sellwoodll2000l ). As an ex- 
ample, Figure [7] shows the /i4 profile along the major axis 
slit. The square symbols are the target data computed us- 
ing higher order Jeans equations without seeing, the circles 
are computed the same way but including PSF effects. The 
lines are the temporally smoothed profiles of the particle 
model using Npp = 5 to represent the PSF. One sees that 
the heavily seeing-affected central profile is well recovered 
by the model. 

In the dynamical modeling of NGC 4697 we include 
seeing only for our new kinematic data. We represent the 
PSF by a single Gaussian with FWHM = 1.25". For the 
BDI data we do not know the PSF but since the slit cells 
are relatively large, seeing is likely to be negligible. 

3.4 The merit function 

By fitting the particle model to the observables, the weights 
Wi are gradually changed such that the merit function 

1 2 



F = fiS 



(12) 



is maximized, where S is a profit function and measures 
the quality of the fit. The parameter fi controls the relative 
contribution of the profit function to F; incrementing ^ in- 
creases the infiuence of S in equation (|12|l . The x^ statistics 
is computed as usual 



(13) 



where Aj — {yj — Yj) /a{Yj). yj is a model observable {e.g. 
o.im,k with j = {lm,k}), Yj is the corresponding target and 
©■(Yj) its error. 

For the profit function S, we use the entropy 
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Figure 7. Seeing convolution test, comparing the radial run of 
/i4 along a 2" slit for a spherical target model and its particle 
representation. The squares (circles) were computed for the target 
from higher order Jeans equations without (with) seeing. The 
lines correspond to the particle model including seeing, for which 



the PSF was represented using A'p 



5. The full and dashed 



lines refer to the major axis slit data at positive resp. negative 
radii with respect to the origin. The heavily seeing-affected central 
profile is well recovered by the model. 



Wi \n{wi/wi) 



(14) 



where {wi} are a predetermined set of weights, the so-called 
priors (here equal for all particles) . The entropy term pushes 
the particle weights to remain close to their priors (more 
specifically, close to Wi/e). This implies that models with 
large ^ will have smoother distribution functions than those 
with small ^. The best choice for fj, depends on the obser- 
vational data to be modeled, e.g. their spatial coverage, on 
the phase-space structure of the galaxy under consideration, 
but also on the initial conditions. For the dataset at hand, 
the best value of fj, will be determined in Section [4.11 



3.5 Discrete PNe velocities 

The likelihood of a model fit to photometric as well absorp- 
tion line kinematic data is measured by the standard 
statistics given in equation (|13|) . To treat discrete PN ve- 
locity measurements the same way, we must bin them to 
estimate the underlying mean v and o fields. This gives 
the corresponding model observables as discussed in section 
[3X2I 

As an alternative, one can measure the likelihood of a 
sample of discrete velocities Vj and positions Rj = {xj,yj) 
on the sky via 



£ = ^ln£j 



(15) 
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using the likelihood functio n for a single PN 
iRomanowskv fc Kochanekll200ll ) 



(16) 



where aj is the error in velocity and dL/dv^ is the LOSVD 
assuming as before that the line-of-sight is along the z-axis. 

We can then add equation (|15p to the function F given 
in equation (I12|l and maximize 



F+ ^F + C 



(17) 



with respect to the particle weights Wi. Hence, we obtain 
an additional contribut ion to the force-of-change as given in 
Ide Lorenzi et al] (|2007l ). We will now derive this extra term. 
Let us consider the selection function 



5ji = 



1 if {xi, yi) e Cj 
otherwise. 



which assigns particle weights to a spatial cell Cj, which 
contains the j-th PNe. We can then write dL/dvz at position 
j as 



with 



(18) 



(19) 



and 5{x) being the standard delta function. Hence, equation 
l|16p can be expressed in terms of the particles via 



' h 
with 



Finally, we find for the additional term in the FOC 



dm 



2tt 



(20) 



(21) 



(22) 



where the sum runs over all individual PNe. For small errors, 
the dwi/dt from the likelihood term is positive for particles 
with Vj = Vz,i, but reduces the weights of the other par- 
ticles and hence drives the LOSVD to peak at Vj. In the 
implementation, we replace and jCj with the correspond- 
ing temporally smoothed quantities. 

When we use this method to account for the PN veloc- 
ities in NGC 4697, we adopt the grid defined in Figure [6] 
by the dotted lines, including the innermost and outermost 
full ellipses. In this way, we assign each of the 762 PNe to a 
cell Cj. It follows, that more than one 
spatial bin, but this is not a problem. 



3.6 Efficient mass-to-light estimate 

It is common practice to evolve N-particle systems in in- 
ternal units (lU), in which the gravitational constant and 
the units of length and mass are set to unity, and to scale 
the system to physical units (PU) a posteriori to compare 



with galaxy observations. Similarly, the velocities of a sys- 
tem with mass-to-light ratio T of unity may be scaled to any 
T via vpu = 7 vm where 7 oc \/T and vp\j and «iu are the 
velocities in physical and internal units, respectively. It fol- 
lows that the kinematic observables of the model and hence 
also can be regarded as functions of T. Equation p3|l 
then reads 



X 



(23) 



In the following we will only consider the luminosity- 
weighted Gauss-Hermite moments as given in equation ((Sjl 
and neglect the PNe kinematic constraints. Taking the par- 
tial derivative with respect to T of equation (|23|) leads to 



19x! 

2 dr 



E Aj(T) db„,p 

j 



j = {n,p} 



(24) 



where -B„,p is the target observable and a{Bn,p) its error. We 
define a force-of-change (FOC) for the mass-to-light ratio T 



dT ^_ V 

dt or 

which equals 

f ^-.T^2A.(T)-^^^m 



with 
dT 



Tapa{Bn 



dT 



du„{x) 



dx 



(25) 



(26) 



(27) 



where we used dv^^i/dT = Vz,i/2T for Vz,i given in physical 
units, and j = {n,p}. The line-of-sight is along the z axis. 
In practice, we use the temporally smoothed quantities to 
compute the FOC for the mass-to-light ratio. 

In principle, the proposed scheme can be understood 
as a gradient search along the X^C*") curve when simultane- 
ously the particle model is fitted to the observational con- 
straints. Hence the same NMAGIC run allows us to estimate 
T as well. We test the scheme and illustrate its accuracy in 
Section [l] 



3.7 Initial conditions 

As initial conditions for NMAGI C, we g enerate a par ticle re- 
alizat ion of a spherical 7-model (jPehnen 1993; Caroll o~et al.l 
Il995h mad e from a distribution f uncti on (DF) using the 
method of lOebattista fc SellwoodI (|200a ). The model con- 
sists of A'^ = 5 X 10^ particles and has 7 = 1.5, scale length 
a — 1 and rmax = 40. When scaled to NGC 4697 one unit 
of length corresponds to 2.3810 kpc, i.e., this model has 
Re = 3.8 kpc. 

In some cases, we have found it useful to give the ini- 
tial particle system some angular momentum about an axis 
of symmetry. For axisymmetric stellar systems, the den- 
sity is determined through the even part in of the DF 
(|Lvnden-Belllll962l ). Thus the component of the angular 
momentum of a particle along the symmetry axis may be 
reversed without affecting the equilibrium of the system. 
iKalnaid |l973) showed, however, that a discontinuity at 
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Lz = can affect the stability of the particle model. There- 
fore, if desired, we switch retrograde particles with a prob- 
ability 


5 





1 1 1 

major axis ! 


^po^2^'^2. (28) 
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which ensures a smooth DF. 
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4 TESTING THE MODELLING WITH 
ISOTROPIC ROTATOR TARGETS 
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In this section, we use axisymmetric, isotropic rotator mod- 
els with known intrinsic properties to determine the optimal 
value of the entropy "smoothing" parameter p in equation 
(|12p. and to test our procedure for determining the optimal 
mass-to-light ratio simultaneously with modelling the data. 
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4.1 Entropy parameter 

Our approach to de t ermin e suitable values for /j is similar 
as in iGerhard et al] (| 19981 ) : Ixhomas et all l|2005h . We first 
generate a "mock" kinematic data set from an isotropic rota- 
tor model whose information content (number and density 
of points, errors) is similar as for the real data set to be 
modelled. To this data set we perform a sequence of particle 
model fits for various ^, and determine the values of ji for 
which (i) a good fit is obtained, and (ii) the known intrin- 
sic velocity moments of the input "mock" system are well 
reproduced by the corresponding moments of the final parti- 
cle model. Using an isotropic rotator model for this purpose 
here makes sense, b ecause such a mode l is a fair representa- 
tion of NGC 4697 jSinnev et al.lll990l '). We have chosen to 
describe the luminosity density of the mock galaxy by one 
of the flattened 7- models of De hnen fc Gerhard (,199 3) , 



j(m) = 



(3-7)i 



47rg mT(m-fa)4-7' 



(29) 



Here L and a are the total luminosity and scale radius, w? = 
R' + {z /q) , and q is the fiattening. The parameters are 
chosen such that the surface brightness closely resembles 
that of NGC 4697, i.e., q = 0.7, 7 = 1.5, L = 2 x 10^° Lq,r 
and a = 2.5 kpc, which corresponds to a ~ 49" at a distance 
of 10.5 Mpc. Figure [8] shows a comparison of the surface 
brightness of NGC 4697 with the mock galaxy projected 
under i = 80°. The major and minor axis surface brightness 
profiles are well approximated by the 7-model, except for 
some differences at larger radii, so we will use this model for 
the calibration of ^. 

We determine mock kinematic profiles from internal ve- 
locity moments, obtained fr om higher-order Jeans equations 
l|Magorrian fc Binnevl[l994l ') in the self-consistent potential 
generated by the density of equation (|29[) for a mass-to-light 
ratio T = 5. Before calculating the line-of-sight velocity pro- 
files, the velocity moments are slit-averaged to account for 
the observational setup of the kinematic slits given in Sec- 
tion [23] We add Gaussian random variates to the isotropic 
rotator kinematics with la dispersion corresponding to the 
respective measurement error in NGC 4697 at that position. 
Figure [5] shows a comparison of our new kinematic data for 
NGC 4697 with the isotropic rotator mock data, along the 
galaxy's major axis. 
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Figure 8. Comparison of the surface brightness of NGC 4697 
(solid lines) with the 7-model described in the text and seen under 
i = 80° inclination (dotted lines). Top: Surface brightness profile 
along the major axis, (b) Along the minor axis. 
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Figure 9. Comparison of v, a, /13 and /14 of NGC 4697 and an 
isotropic rotator model with approximately the same projected 
surface brightness as the galaxy. The filled circles show our new 
kinematic data for NGC 4697 from Section 12.31 the star sym- 
bols show the isotropic rotator mock data, and the dashed red 
lines show the underlying smooth model kinematics, all along the 
major axis. 



We do not construct mock PNe data for inclusion in 
the entropy tests, but we need to compute the photometric 
observables to construct a complete observational data set. 
We expand the luminosity distribution of equation (|29|l in a 
spherical harmonics series (c/. Section [3. 2. l|l on a radial grid 
with 40 shells at radii r^. The radii are quasi- logarithmically 
spaced with rmin = 1.0" and rmax = 700". We use the lumi- 
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nosity on radial shells Lk = \/47ryloo,fc and the higher order 
coefficients A2o,k, ^22, k, Aae.k and Aso,k to constrain 
the luminosity distribution of the particle model. The m 7^ 
terms are set to zero to force the models to remain nearly 
axisymmetric, i.e., within the limits set by the Ai m errors. 
We assume Poisson errors for the Lfci a{Lu) = LkL/N 
where N is the total number of particles used in the particle 
model and L is the total luminosity of the system. To esti- 
mate the errors in the higher order luminosity moments, we 
use Monte-Carlo experiments in which we generate particle 
realizations of a spherical approximation of the density field 
of the target system with 5 x 10^ particles, which is the same 
number as in the x^M2M models for NGC 4697. 

We then construct self-consistent particle models for 
the isotropic rotator target in a three step process, us- 
ing the mock data as constraints, (i) Density fit: we start 
with the spherical initial conditions described in Section 
13.71 and evolve them using NMAGIC to generate a self- 
consistent particle realization with the desired luminosity 
distribution (7-particle model), fitting only the luminosity 
constraints, (ii) Kinematic fit: because the target galaxy has 
a fair amount of rotation, it is worth starting the kine- 
matic fit from a rotating model. Hence, following Section 
13.71 we switch a fraction of retrograde particles in the 7- 
particle model to prograde orbits, using po = 0.3 and 
Lo — I/circ(0.03-Re). This rotating system we then use as a 
starting point to construct a series of self-consistent dynami- 
cal isotropic rotator models, by fitting the target photometry 
and kinematics for different values of /i. For each model, we 
evolve the particle system for ~ 10^ NMAGIC correction 
steps while fitting the complete set of constraints. During 
this correction phase, the potential generated by the parti- 
cles is updated after each correction step, (iii) Free evolution: 
to ensure that any correlations which might have been gen- 
erated during the correction phase are phase-mixed away, we 
now keep the potential constant and evolve the system freely 
for another 5000 steps, without further correction steps. For 
reference, 5000 of these steps correspond to ~ 20 circular 
rotation periods ("dynamical times" )at in spherical ap- 
proximation. 

The results are presented in Figure [10] The lower panel 
shows the quality of the fit as a function of fi, both in terms 
of normalized values and in terms of the merit function 
F from equation (|12|l . The upper panel shows the rms A 
relative difference between the internal velocity moments of 
the isotropic rotator input model and those of the parti- 
cle models reconstructed from the mock kinematics. For the 
particle models intrinsic velocity moments are computed by 
binning the particles in spherical polar coordinates, using 
a quasi- logarithmic grid with 20 radial shells bounded by 
»'min = 0.01" and rmax = 200", 12 bins in azimuthal angle 
<j), and 21 bins equally spaced in sin 6. The rms A shown 
in Figure [10] is obtained by averaging over all grid points 
in the radial region constrained by the data {R 5C 1.5-Re). 
The minimum in rms A determines the value of ^ for which 
the model best recovers the internal moments of the input 
model. This occurs at ^ ~ 10^, and the value of rms A at 
the minimum is ~ 1.4% . For larger (smaller) ^, the rms A 
is larger because of oversmoothing (excess fluctuations) in 
the model. 

/N values are given in the lower panel of Fig. [10] 
for all (photometric and kinematic) data points, and for 
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Figure 10. Entropy tests. Top: Deviation rms A{fi) of the parti- 
cle models from the isotropic rotator internal velocity moments. 
The point for the rightmost value of fi is at a large value of A 
outside the diagram. Bottom: deviation per data point of the 
particle model fit to the photometric and kinematic target ob- 
servables (open circles) and to the kinematic observables alone 
(filled squares), as a function of entropy parameter /i. The trian- 
gles show the same dependence for the merit function {—F), cf. 
equation H12| l. The starred symbol indicates the value of /i chosen 
for the subsequent modelling. 



the kinematic data points alone. Generally the /N for 
the photometric points is significantly better than for the 
kinematic points, because (i) the Aim come from averages 
over many particles, thus have little noise, and we have not 
added Gaussian variances, and (ii) all particles contribute 
to the Aim force-of-change at all timesteps, so the Aim are 
weighted strongly during the evolution. The kinematic 
per data point in the lower panel is of order 0.5 for a large 
range of ^ and then increases starting from ^ > 300 to 1 at 
^ ~ 5 X 10^, whereas —F already increases around ^ > 100. 

Some results for the isotropic rotator dynamical models 
obtained with — lO'^, = 10^, n = 5 x 10^ are presented 
in Figures [11] and 1121 Figure [11] shows a comparison of the 
target kinematics with the kinematics of the self-consistent 
particle models along the major axis slit. Note the excellent 
fit of the central velocity gradient and velocity dispersion 
dip, for all values. However, the models with higher ^ 
begin to fail matching the target data at the largest radii. 
This is because the number of data points decreases with 
radius, whereas the number of particles and hence entropy 
constraints is roughly proportional to luminosity Lk, i.e. 
changes much more slowly with radius. The result is that 
the constraints from the data become relatively weaker at 
larger radii. The entropy term tries to enforce a dynamical 
structure related to the initial particle model, in which all 
particles have equal weights. In the present case this works 
to cause a bias against both fast rotation and anisotropy. 
This first becomes apparent where the relative statistical 
power of the data is weakest, i.e., at large radii. 

Because our goal in this paper is to determine the range 
of potentials in which we can find valid dynamical models for 
NGC 4697, we need to ensure that the answer to this ques- 
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Figure 11. Particle model fits to the isotropic rotator mock kine- 
matic data along the model's major axis. The points with error 
bars show the target data and the lines represent the model kine- 
matics. The model data points are averages over the same slit cells 
as the target data (see Fig.[5ll, and are connected by straight line 
segments. The model v, a are determined via eqs. (|9]l and HIOI I 
and /i3, h4 are the fitted values based on the observed scale pa- 
rameters and (Tp. The full, dotted and dashed lines correspond 
to the models obtained with = 10'^, ii = 10'', and ^ = 5 X 10'^, 
respectively. The red dashed line shows the fi = 100 model kine- 
matics 20 dynamical times later, reflected with respect to the 
origin, and obtained from direct fitting of the model line pro- 
files. This proves that this model is accurately axisymmetric and 
stationary; see text for details. 



tion is not biased by overly strong entropy smoothing in the 
galaxy's outer regions. Thus in the modelling in Section [5] 
we will conservatively choose /x = 100 for the smoothing pa- 
rameter (indicated by the starred symbol in Fig. llOp . Sim- 
ilar caution is common practice in determinin g black hole 
masses in galaxies (e.g.. ICebhardt et al.ll2003l ). The result- 
ing dynamical models will then be somewhat less smooth 
than could be achieved, but this price is rather modest; be- 
tween fi — 10^ and its minimum value at ^ = 10^, the rms A 
in Fig. [TOl decreases from ~ 1.7% to ~ 1.4%, i.e., by ~ 15%. 
Certainly it would not be appropriate to rule out potentials 
in which the solutions differ by this degree in smoothness. 

Using /i — 100 in the modelling leads to a slight over- 
fitting of the slit kinematic data, especially for the higher 
order kinematic moments which themselves take only val- 
ues of order percents. It is worth pointing out that, con- 
trary to first appearances from Fig. 1111 this implies neither 
that these models are not axisymmetric, not that they are 
out of equilibrium. The model kinematics shown in Fig. [11] 
are obtained after 20 dynamical times of free evolution in 
the axisymmetric potential, so are thoroughly phase-mixed. 
The model data points shown are averages over the same 
slit cells as the target data (see Fig. [5]l , and are connected 
by straight line segments. The plotted v, a are determined 
via eqs. (|9]) and (|10p and /13, are the fitted values based 
on the observed scale parameters Vp and CTp. The red dashed 
line in Fig. [TT] also corresponds to the fi — 100 model, but 



has been determined as follows; (i) from the particle dis- 
tribution after a total of 40 dynamical times of free evo- 
lution; (ii) from a mirror-symmetric set of slit cells, with 
respect to the major-axis slit shown in Fig. [5] (iii) using 
the {v, a, h^, hi) parameters obtained by direct fits to the 
model line profiles, and (iv) finally refiecting the kinematics 
so obtained anti-symmetrically with respect to the origin. 
The excellent agreement between this curve and the orig- 
inal major axis kinematics of this model in Fig. [TT] shows 
that (i) the /i = 100 model is a true equilibrium, (ii) it is 
accurately axisymmetric, and (iii) the left-right differences 
in the kinematics in Fig. [TT] are due to slightly different 
slit cell averages over the model on both sides. That these 
averages can be slightly different is made possible by low- 
level (axi-symmetric) structure in the model consistent with 
the slight under-smoothing for this value of fi. What hap- 
pens is that the algorithm adds a few near-circular orbits in 
the relevant radial ranges. When added to the correspond- 
ing model LOSVDs and averaged over the asymmmetric slit 
cells these orbits change the kinematic moments hn, n ^ 3 
at the ~ 0.01 level so as to improve the agreement with the 
observed major axis hs, /14. In the other slits the models 
interpolate more smoothly between points when needed be- 
cause fluctuations in the particle distribution to follow local 
kinematic features are less easily arranged; see the corre- 
sponding figure for NGC 4697 in Section [5] 

A comparison of the internal velocity moments of the 
input model and the particle model in the equatorial plane 
is presented in Figure [12] The figure shows an, and cr^, 
followed by v^. The last panel displays the anisotropy pa- 
rameter Pe = 1 — (Jg /cTr , which is zero for the input isotropic 
rotator model. Within the radial extent of the target data, 
the internal moments of the input model are well repro- 
duced; outside this region, where the model is poorly con- 
strained by the input data, the particle model increasingly 
deviates from the target. Indeed, if we add PN velocity data 
in this test, the corresponding particle model agrees with 
the internal moments of the input model out to larger radii. 



4.2 Mass-to-light ratio 

We will now use such isotropic rotator models to explore 
how accurately we are able to recover the input mass-to- 
light ratio, given the spatial coverage of the data. Further, 
we will test our new procedure, described in Section [3.61 for 
estimating the mass-to-light ratio efficiently. As input mod- 
els we take both the self-consistent isotropic rotator model 
described above and a model constructed in the same way 
but including a dark matter halo. The halo potential is of 
the form of equation ([TJ, with r\) = 190 arcsec (9.7kpc), 
vo — 220 km s~^ and = 1.0. The mass-to-light ratio of 
the stars in both input models is fixed to T = 5. 

The results for a "classical" approach, in which we fit 
a dynamical particle model to the data for different values 
of T, are presented in Figure 1131 which shows the quality 
of the fit as function of T for the self-consistent case. The 
input value of T is recovered well. The results with the new 
procedure presented in Section [3^ are summarized in Figure 
1141 The figure shows the evolution of the mass-to-light ratio 
as a function of time during NMAGIC model fits. Models 
for both the self-consistent input galaxy and for the target 
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Figure 12. Comparison of tlie intrinsic velocity moments in the 
equatorial plane of the axisymmetric isotropic rotator particle 
model and target model. The points represent the target system 
and the lines correspond to the final particle model for fi = 100, 
averaged over azimuth. The dashed vertical lines show the max- 
imum radial extent of the minor axis (left line) and major axis 
target kinematic data (right line). At larger radii the particle 
model is poorly constrained by the input target data. 
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Figure 13. Quality of the particle model fit to the self-consistent 
isotropic rotator input model, as a function of assumed mass-to- 
light ratio T. values per data point are given for the parti- 
cle model fit to the photometric and kinematic target observ- 
ables (open circles) and to the kinematic observables alone (filled 
squares). The triangles correspond to the measured merit func- 
tion F. The input mass-to-light ratio is T = 5. All models are 
built using fi = 100. 
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isotropic rotator in a potential in- 
lines). The input mass- to-light ratio 
given in terms of elapsed time steps 
to ~ 40 circular rotation periods at 



model including a dark halo potential are shown, with both 
low and high initial choices of T. The tests show that for the 
self-consistent case the input mass-to-light ratio is recovered 
very well. The uncertainties are slightly larger when a dark 
matter halo is included, but the maximum fractional error 
is less than three percent. We conclude that the new scheme 
works very well and that T is recovered within a few percent 
(for the amount and quality of data used in the present 
work). The advantage of the new method is its efficiency, 
only one run is needed to estimate T instead of order 10, 
but at the cost of not knowing the shape of a-s a function 
of T near the minimum, i.e., the confidence interval. 



5 DYNAMICAL MODELS OF NGC 4697 

After these tests we are now ready to use NMAGIC for 
constructing axially symmetric dynamical models of NGC 
4697. We investigate self-consistent models as well as mod- 
els including dark matter halos, and fit the photometry, slit 
kinematics and PNe data. Our aim in this paper is not to 
attempt to constrain detailed halo mass profiles, but only to 
ascertain whether a dark matter component is allowed, or re- 
quired, by the kinematic data. Thus we investigate a simple 
sequence of potentials A to K which include the contribution 
from the stellar component and a halo potential as in equa- 
tion IT]), with parameters given in Table [1] The parameters 
are chosen to result in a sequence of circular speed curves 
ranging from falling according to the distribution of stars 
to nearly flat over the whole range of radii. This sequence 
is shown in Figure [15] all these circular velocity curves are 
computed in the galaxy's equatorial plane and include the 
stellar component with the respective best-fitting mass-to- 
light ratio as given in Table [1] 
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vq /km s 


Q<t> 










-F 


T 


A 








1.0 


0.453 


0.0323 


0.900 


0.968 


370.2 


5.78 


B 


5.76 


80 


1.0 


0.415 


0.0254 


0.828 


0.884 


343.9 


5.74 


C 


5.76 


120 


1.0 


0.439 


0.0343 


0.877 


0.784 


358.6 


5.71 


D 


4.32 


160 


1.0 


0.404 


0.0288 


0.816 


0.610 


333.7 


5.58 


E 


4.32 


190 


1.0 


0.404 


0.0244 


0.826 


0.520 


332.8 


5.49 


F 


4.32 


210 


1.0 


0.386 


0.0229 


0.791 


0.476 


320.0 


5.45 


G 


4.32 


210 


0.8 


0.382 


0.0203 


0.785 


0.439 


315.4 


5.46 


H 


2.88 


210 


0.8 


0.376 


0.0232 


0.773 


0.397 


310.2 


5.28 


J 


4.32 


250 


0.8 


0.383 


0.0242 


0.786 


0.377 


313.7 


5.34 


K 


2.88 


250 


0.8 


0.377 


0.0212 


0.771 


0.506 


309.6 


5.10 



Table 1. Table of model parameters and fit results. Columns (l)-(4) give the model code and the parameters ro, vo and used in 
equation JTJl for the respective dark halo potential in this model. The next four columns list the values per data point, for all observables 
[column (5)], and for the density constraints, slit kinematic observables, and PN observables (data set PNDl) separately [columns (6)-(8)]. 
Column (9) gives the numerical value of the merit function in equation 1121 and column (10) the final (r-band) mass-to-light ratio. The 
respective number of constraints are A'^ = 1316, Naim = 680, N^i = 604, Npi^ = 32. 
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120 
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0.419 


0.0331 


0.853 
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5.72 
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160 


1.0 


0.406 


0.0314 


0.828 


2028.3 


2357.9 


5.60 
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4.32 


190 


1.0 


0.391 


0.0271 


0.801 


2026.3 


2344.9 


5.54 


F 


4.32 


210 


1.0 


0.402 


0.0304 


0.820 


2025.6 


2350.1 


5.49 


G 


4.32 


210 


0.8 


0.396 


0.0232 


0.815 


2024.8 


2343.9 


5.48 


H 


2.88 


210 


0.8 


0.373 


0.0245 


0.766 


2026.3 


2329.2 


5.31 


J 


4.32 


250 


0.8 


0.374 


0.0203 


0.773 


2025.6 


2329.0 


5.37 


K 


2.88 


250 


0.8 


0.369 


0.0198 


0.763 


2030.8 


2329.9 


5.14 



Table 2. Table of model parameters and fit results, similar to Table [T] but with all models computed using the likelihood scheme for the 
PNe as discrete kinematic tracers. Columns (8) and (9) now give the likelihood of the PN data set C and the merit function including 
C [equation [T7] . The other columns are equivalent to those in Table [l] 



To construct the models, we proceed as in Section U 
First, we compute the photometric observables. We expand 
the deprojected luminosity distribution of NGC 4697 in 
a spherical harmonics series on a grid of 40 shells in ra- 
dius, quasi-logarithmically spaced with rmin ~ 1-0" and 
Tmax = 700" . As obsetvablcs we use the luminosity on radial 
shells Lk and the higher order coefficients A2o,k, ^22, fe, ■ ■ ■, 
^66, fc and Ago.fc, at radii r^. The m 7^ terms are set to 
zero to force the models to remain nearly axisymmetric, i.e., 
within the limits set by the specified Aim errors. Because 
the photometry is not seeing-deconvolved, for the innermost 
two points {R < 3" ) we only fit the ^00 term. Errors for the 
luminosity terms Aim are estimated by Monte Carlo simula- 
tions as in Section [4.11 As kinematic constraints we use the 
luminosity weighted Gauss-Hermite moments from the slit 
data, and the PNe kinematics, either represented by binned 
line-of-sight velocity and velocity dispersion points, or as 
discrete velocity measurements; see Sections 12.3.21 and 13.51 

Again we fit particle models in a three step process. 

(i) First, we start with the spherical particle model de- 
scribed in Section|3T7]and evolve it using NMAGIC to gener- 
ate a self-consistent particle realization with the luminosity 
distribution given by the deprojection of the photometry. 

(ii) Because NGC 4697 shows significant rotation, we then 
switch retrograde particles similarly as in Section 14.11 us- 



ing po = 0.3 and Lo — Lcirc(0.03i?e). The resulting rotating 
particle model (hereafter, model RIC) is used as a starting 
point to construct a series of dynamical models by fitting 
the photometry and kinematics in different halo potentials, 
as follows. For every dark matter halo from Table [T] we first 
relax RIC for 5000 time steps in the total gravitational po- 
tential, assuming a mass-to-light ratio of 5.74. For reference, 
10000 time steps in the self-consistent potential correspond 
to ~ 40 circular rotation periods at Re in spherical approx- 
imation. After this relaxation phase, we evolve the particle 
system for ~ 10^ NMAGIC correction steps while fitting the 
complete set of constraints. During this correction phase, the 
potential generated by the particles is updated after each 
correction step but the dark matter potential (if present) 
is constant in time, (iii) Subsequently, we keep the global 
potential constant and evolve the system freely for another 
5000 steps, without further correction steps. Models A, D, 
G and K were in addition evolved for a further 10000 steps 
with all potential terms active, to confirm that the modest 
radial anisotropy required in these models does not lead to 
dynamical instabilities. 

To make sure that the results are not biased by the way 
we incorporate the PNe data, we have constructed three 
models in most halo potentials. Each time the PNe data are 
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Figure 15. Circular velocity curves of the potentials used in the 
modelling, including the self-consistent model A (dashed line), 
and a sequence of dark matter halos (solid lines). The lines at 
r/Re = 7 run from model A (bottom) to K (top), with models F 
and G represented by the same curve; cf. Table [T] 



Figure 16. Comparison of the photometric constraints with the 
final models A (self-consistent, dashed), D (dotted), and G (full 
line). The points correspond to the target input data. From top 
to bottom: luminosity on radial shells profile L^. = \/47rAoo, and 
normalized A20 and j44o profiles. 



represented differently, using the binnings PNDl, PND2, or 
the likelihood method. 

The quality of the fit for different halo models can be 
characterized by the quantity F defined in equation (|12p 
or (|17p and is given in Tables [T] and [2] In addition, the 
value of P^r data point is also shown, globally and for 
each data set separately, as are the stellar mass-to-light ra- 
tios. For the same reasons as for the isotropic rotator test 
models, the density constraints are very accurately fit. The 
slit kinematics are typically fit within about 0.9a per point, 
slightly better than required. This is due to the relatively low 
value used for the entropy smoothing, needed not to bias the 
range of allowed potentials by the imposed smoothing. The 
PN values indicate that the PN data are consistent with 
all models. 

The likelihood values for the models reported in Table 
[2] would formally allow us to exclude a large fraction of the 
halo potentials tested, using AL = 0.5. However, these like- 
hhood differences depend to a significant degree on very few 
PNe in the wings of the LOSVDs. To assess this we have 
compared histograms of PN velocities with the LOSVDs of 
these models in cones along the major and minor axes. A 
test shows that all model LOSVDs except for the model 
without dark halo are formally consistent with the PNe ve- 
locity histograms and associated Poisson errors. However, 
the low dark matter models show systematic deviations in 
form from the more fiat-topped data histograms, which we 
judge significant for models A to C. Also, for these models 
the LOSVD a-nd the likelihood are correlated, which is 
not the case for the more massive halo models. On the basis 
of these findings we believe the likelihood results can prob- 
ably be trusted for ruling out models A-C, but not for any 
of the more massive halo models. Thus we estimate that 
the range of circular velocities at 5Re consistent with the 
data is Vc{5Re) ^ 200km s~^, with the best models having 



Vc{5Re) — 250km s~^. To reach more stringent constraints 
will require PN velocities at even larger radii. 

Figures 1161 1171 1181 present results from some of these 
models, comparing the stars-only model A and the three 
halo models D, G, and K to the data. Figure [TS] shows the 
comparison of models A, D, G with the photometric con- 
straints. The model lines match the target data points per- 
fectly, in accordance with the very small Xaim/^aim values 
in Table[T] Figure [l7]compares the projected absorption line 
kinematics of the three models with our measurements and 
the BDI data. The fits are generally excellent. Along the ma- 
jor and minor axes one can see how the models have found 
compromises to deal with asymmetries of the data on both 
sides of the galaxy, and slight discrepancies between our and 
the BDI data, e.g., in the region around ±10" along the ma- 
jor axis. As for the isotropic rotator, the major axis higher 
order moments in Fig. [T7] are even somewhat overfitted; see 
the discussion in Section [4.11 

Figure [THI compares the final A, D, G, and K models 
with the PNe kinematic constraints along the major axis 
(left) and minor axis (right); on each axis we show mean 
velocity (top) and velocity dispersion (bottom). The model 
curves in Fig. [T5] and the x^ per data point values in Table 
[T] are computed for PN dataset PNDl. There was no dif- 
ference between these values and those obtained with PN 
dataset PND2 in all cases where we modelled both. The two 
additional lines in the panels of Fig. [TS] show the mean ve- 
locities and velocity dispersions for the variants of models A, 
K obtained with the likelihood scheme for the PNe (see Sec- 
tion |33J, computed by binning the particles in these models 
a posteriori in the same bins as for dataset PNDl. While 
there is little difference for model K, the likelihood variant 
of model A fits the observed PN data points actually better 
than the original model A based on the PNDl data. 

Overall, this figure illustrates that with increasing halo 
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Figure 17. Comparison of models A, D, and G to tlie absorption line kinematic data along the major axis (top left), minor axis (top 
right), the slits parallel to the major axis (bottom left), and the slits parallel to the minor axis (bottom right). Full and starred data 
points show our new data and the BDI data, respectively. The model data points are averages over the same slit cells as the target data 
(see Fig. (Sj, and are connected by straight line segments. Linestyles for the models are the same as in Fig. 1161 



mass the fit to the PNe kinematic data improves slightly. 
Models G and K bracket the best-fitting models to the 
binned data in Fig. 1181 However, also model A without dark 
matter still has a Xpn / ^pn .just below one, despite being 
systematically a little low in the minor axis dispersion plot. 
When we force the self-consistent model to improve the PN 
data fit at the expense of the slit kinematic data fit, the 
model starts to develop specific anisotropy features at the 
radii of the outer PN data. This suggests that with PN data 
extending to somewhat larger radii, « 400 arcsec [R ~ 6i?e), 
the model without dark halo might start to fail also for the 
binned data in Fig. 1181 consistent with the likelihood and 
LOSVD results. In conclusion, a variety of dark halos with 



Vc{5Re) ^ 200km s ^ are consistent with all the kinematic 
data currently available for NGC 4697. 

Finally, Figure [19] shows the internal kinematics of the 
particle models A, D, G, and K. The upper panels give 
(Tij, 0"^ and (Tz, followed by v^. The last panel displays the 
anisotropy parameter /3e = 1 — ag/a'^, which is zero for an 
isotropic rotator model. All quantities are given as aver- 
ages over the models' equatorial plane. The more massive 
halo models become more radially anisotropic in the outer 
parts in terms of crjj vs. a^, but Pe does not increase beyond 
model D because and increase in parallel while 
remains constant. Thus the additional kinetic energy that 
stars at large radii must have in these models, is hidden in 



Dark matter content and internal dynamics of NGC 4697 17 



major axis 



minor axis 



6 SUMMARY AND CONCLUSIONS 




100 200 300 
R (arcsec) 



50 100 150 200 
R (arcsec) 



Figure 18. Comparison of the PNe velocity and velocity disper- 
sion data (PNDl, points) with models A, D, G, and K. Top left: 
V along the positive major axis. Top right: The same for the mi- 
nor axis. Bottom left: cr along the positive major axis. Bottom 
right: The same but for the minor axis. Dashed, dotted, full, and 
upper dashed lines show models A, D, G, and K; the two dash- 
dotted lines show the variants of models A and K obtained with 
the likelihood scheme for the PNe. 
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Figure 19. Internal velocity moments in the equatorial plane for 
models A, D, G, K (dashed, dotted, full, and dash-triple dotted 
lines, respectively). The vertical dashed lines indicate the radial 
extent of the minor axis slit data, major axis slit data, and PN 
data, from left to right. 



the plane of the sky. Conversely, at small radii the velocity 
dispersions in models G-K are slightly lower, compensating 
for the larger radial velocities of halo stars along the line-of- 
sight to the center. These models have /3 ~ 0.3 at the center, 
which increases with radius and reaches /3 ~ 0.5 at > 2Re. 



In this paper, we have presented new surface brightness mea- 
surements and long slit spectroscopic data for the E4 galaxy 
NGC 4697, and combined these data with existing long slit 
kinematics and discrete PNe position and velocity measure- 
ments to construct dynamical models for this galaxy. The 
combined data set runs from the center of the galaxy to 
about 4.5 effective radii. 

For the first time, we have modelled such a dataset with 
the new and flexible x^-made-to-measure {x'^M.2M) particle 
code NMAGIC. We have extended NMAGIC to include see- 
ing effects and have implemented an efficient method to es- 
timate the mass-to-light ratio T. Tests of this scheme using 
isotropic rotator input models have shown that the method 
recovers T within a few percent both for self-consistent and 
dark matter dominated target galaxies. In addition, we have 
implemented a likelihood scheme which allows us to treat the 
PNe as discrete velocity measurements, so that no binning in 
velocities is needed. The modelling presented in this paper 
shows that the x^M2M/NMAGIC particle method is now 
competitive with the familiar Schwarzschild method. In fact, 
it has already gone further in that the gravitational potential 
of the stars has been allowed to vary in the modelling, the 
mass-to-light ratio has been adapted on t he fly, the stabil- 
ity of the models has been checked, and, in lde Lorenzi et akl 
|2003), NMAGIC has been used to construct triaxial and 
rotating triaxial models. 

Even though NMAGIC does not require any symme- 
try assumptions for the modelling, we have in this paper 
forced the method to generate axisymmetric particle mod- 
els for NGC 4697. Both self-consistent models without dark 
matter, and models following a sequence of circular speed 
curves with increasing dark halo contributions have been 
investigated. The PN data have been used both binned on 
two different spatial grids, as well as with the new likelihood 
scheme, to make sure that the results are not biased by the 
way the PNe data are incorporated. 

Our main astronomical result is that models with a va- 
riety of dark matter halos are consistent with all the data 
for NGC 4697, as long as the circular velocity Vc{5Re) ^ 
200km s~^ at 5Re- These models include potentials with 
sufficiently massive halos to generate nearly flat circular 
rotation curves. These models fit all kinematic data with 
X^/N < 1, including the PN LOSVDs. Models with no dark 
matter are not consistent with the PN LOSVDs, as judged 
from both the LOSVD histograms and their likelihood val- 
ues. Amongst the acceptable models, the more massive dark 
halo models with Vc{5Re) — 250km s~^ tend to fit the data 
slightly better in the sense of lower ^ for both the slit 

kinematics and the PN data, but these variations are small 
and not yet statistically significant. To further narrow down 
the range of acceptable dark matter models would require 
PN velocities at even larger radii than currently available, 
out to an estimated ~ 6J?e from the center. 

Our models diff er from earlier studies performed by 
iMendez et all (|200ll ) in the sense that we generate axisym- 
metric models instead of spherical ones and that our models 
are fiexible with regard to anisotropy. The best-fitting mod- 
els are slightly radially anisotropic, with /3 ~ 0.3 at the 
center, increasing to (i 0.5 at > 2 7 ?^. T his is consistent 
with the value given by iDekel et akl (|2005l ) obtained from 
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merger simulations carried out within the ACDM cosmol- 
ogy framework. 
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APPENDIX A: PHOTOMETRIC AND 
KINEMATIC DATA 



Table El Photometric parameters of NGC 4697 
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° Statistical errors not including systomatics duo to photometric calibration and sky subtraction 

^ Error on the center coordinates from the residual rms of the ellipse fit to the isophotes: Err=rms(fit) / \/N with N 128 the number of fitted points of the isophotes. 
" Error of Fourier coefficients defined as Err= y ' ^ a • 
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Table Al. The kinematics of NGC 4697 along the major axis (P.A.=66°). Positive radii are to the north-east. 



R 


V 


dV 


(J 


da 


^3 








(") 


(km/s) 


(km/s) 


(km/s) 


(km/s) 










U.oO 


1 Q 9 

-10. Z 


U.4t 


loU.U 


n A 

U.4 


-U.UU4 


n nn9 

U.UUZ 


O07 
U.UU ( 


no9 

U.UUZ 


u.y ( 


-oy.o 


n 9 
U.z 


1 7c; 
1 / 0.0 


n A 
u.4 


n 097 
U.UZ / 


n nn9 

U.UUZ 


01 /I 
U.U14 


009 

U.UUZ 


1.0 / 


-01.0 


fl 9 

u.z 


1 74 A 


n A 

u.4 


n DAI 

U.U4± 


n 009 

U.UUZ 


017 
U.Ul / 


009 
U.UUZ 


Z.Z f 


7n 1 

- ( U.l 


U.o 


"1 fiO A 
IDVJ.4 


n /I 

U.-l 


u.U4y 


A09 
U.UUZ 


090 

u.uzy 


009 

U.UUZ 


0. ID 


s*? 1 
-00. 1 


u.o 


IO0.4 


n A 
u.4 


U.Uo 


no9 

U.UUZ 


0'?7 
U.UO ( 


009 
U.UUZ 


A Oft 


SQ Q 

-oy.y 


U.o 


lOU.D 


n /I 
u.4 


u.uyo 


n no9 

U.UUZ 


CX'XA 
U.Uo4 


n 009 

U.UUZ 


0.00 


-yo.cs 


n A 

U.4t 


1 "^7 Q 

10 1 .y 


u.o 


U. lUO 


00*^ 

U.UUo 


n 0*^7 

U.Uo ( 


00"^ 
U.UUO 


7 1/1 

( .141: 


Q7 Q 

-y ( .0 


n /I 

U.4I: 


1 ^19 Q 
lOZ.O 


u.o 


n 1 9 
u. iz 


no9 

U.UUZ 


U.UoO 


n no9 

U.UUZ 


Q 99 


Q4 9 

-y4i:.z 


n A 

U.4t 


1 fi'i "i 
100.0 


u.o 


u.uyo 


n09 

U.UUZ 


01 Q 

u.uiy 


00*^ 
U.UUo 


1 9 nn 

iZ.UU 


-yo.D 


n /I 


IDD.O 


u.o 


u.uytj 


n09 

U.UUZ 


01 n 
u.uiu 


no9 

U.UUZ 


iO.D ( 


-yo.o 


U.4 


1 70 ^ 
1 / U.O 


U.D 


U.UoU 


009 
U.UUZ 


OOA 
U.UU4 


00*? 
U.UUO 


on fiQ 


11/11 
- 1 14. 1 


U.D 


1 fin Q 


u.o 


U. lOU 


00*3 

U.UUo 


09 9 
U.UZZ 


00/1 
U.UU4 


97 A9 
Z ( .DZ 


1 HQ 1 

- luy . 1 


U.D 


10 1 .0 


u.o 


fl 1 AO 
U. 14U 


00*^ 
U.UUO 


01 
-U.UID 


OOA 
U.UU4 


QQ n7 
00. U ( 


1 1 1 

-110.1 


n 7 

U. ( 


1 f^l 

101. y 


1 n 
l.U 


fl 1 99 
U. IZZ 


00 /I 
U.UU4 


00*3 
-U.UUo 


U.UUO 


Oo.UO 


1 HQ « 

-lUo.D 


1 n 

l.U 


1 /I Q 9 


1 A 
1.4 


A 1 99 
U. IZZ 


OOfi 
U.UUD 


090 

-U.uzy 


007 
U.UU / 


Q9 c:o 

yz.oz 


111 
-111.0 


9 Q 
z.y 


1 /in Q 
i4u.y 


Q 7 
0. / 


n m 7 

-U.Ul / 


090 

u.uzu 


nc9 

-U.UoZ 


ni 
U.UIO 


n 99 
-u.zz 


7 n 
( .u 


n /I 

U.4I: 


1 7Q « 

1 ( y.o 


u.4 


-U.UZ4 


009 
U.UUZ 


n 01 c; 

U.UIO 


no9 

U.UUZ 


-U.Ol 


Q1 Q 


U.o 


1 77 n 
1 / 1 .u 


n /I 
u.4 


n n/19 

-U.U4Z 


009 
U.UUZ 


n 09K 

U.UZD 


n noi 

U.UUl 


1/11 


C^l Q 

01. y 


U.o 


1 7Q Q 
1 / 0.0 


n /I 
u.4 


-U.UDl 


009 
U.UUZ 


09Q 

U.uzy 


noi 

U.UUl 


-2.11 


70.1 


0.3 


169.0 


0.5 


-0.072 


0.002 


0.045 


0.002 


-3.00 


85.6 


0.3 


163.4 


0.5 


-0.100 


0.002 


0.052 


0.002 


-4.00 


92.8 


0.3 


161.8 


0.3 


-0.112 


0.001 


0.034 


0.002 


-5.19 


93.6 


0.3 


162.7 


0.3 


-0.112 


0.001 


0.041 


0.002 


-6.69 


98.0 


0.3 


160.3 


0.3 


-0.135 


0.001 


0.043 


0.002 


-8.57 


96.0 


0.3 


159.1 


0.3 


-0.114 


0.001 


0.034 


0.002 


-11.05 


98.3 


0.3 


165.8 


0.3 


-0.131 


0.001 


0.030 


0.002 


-14.32 


93.2 


0.3 


172.2 


0.3 


-0.126 


0.001 


-0.014 


0.002 


-18.79 


103.9 


0.3 


172.1 


0.3 


-0.145 


0.001 


-0.011 


0.002 


-24.99 


112.6 


0.4 


168.1 


0.3 


-0.170 


0.002 


0.005 


0.002 


-34 


120.1 


0.5 


160.0 


0.5 


-0.168 


0.002 


-0.026 


0.003 


-48.87 


113.7 


0.7 


152.5 


0.6 


-0.132 


0.004 


-0.018 


0.004 


-76.65 


115.2 


1.9 


153.1 


2.1 


-0.087 


0.012 


-0.004 


0.010 



Table A2. The kinematics of NGC 4697 along the minor axis (P.A.=156°). Positive radii are to the south-east. 
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Figure Al. Isophotal parameters of NGC 4697 as a function of the logarithm of the semi-major axis distance in arcsec. The radial 
profiles of the R-band surface brightness, third, fourth, and sixth cosine Fourier coefficients (03, 04, and ag), and a;— coordinate of the 
center Xcen are plotted in the left panels (from top to bottom). The surface brightness is shown along the major (upper profile) and minor 
axis (lower profile). The radial profiles of the position angle (PA), ellipticity (1 — 6/0), third fourth, and sixth sine Fourier coefficients 
(63, 64, and bo), and y— coordinate of the center {Ycen) are plotted in the right panels (from top to bottom). 



